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Abstract 



On a Ginzburg-Landau model for the macroscopic behaviour of a shape memory aUoy is pro- 

posed. The model is one-dimensional in essence, in that we consider the effect of the marten- 
sitic phase transition in terms of a uniaxial deformation along a fixed direction and we use a 
scalar order parameter whose equilibrium values describe the austenitic phase and the two 
martensitic variants. The model relies on a Ginzburg-Landau free energy defined as a func- 
tion of macroscopically measurable quantities, and accounts for thermal eflFects; couplings 
W between the various relevant physical aspects are established according to thermodynamic 

^ consistency. The theoretical model has been implemented within a finite-element framework 

g and a number of numerical tests are presented which investigate the mechanical behaviour 

I of the model under diflFerent conditions; the results obtained are analysed in relation to 

jH experimental evidences available in literature. In particular, the infiuence of the strain-rate 

O and of the ambient conditions on the response of the model is highlighted. 

""^ Keywords: phase field model, shape memory alloys, strain-rate sensitivity, 

i-H thermo-mechanical coupling 
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VO 1. Introduction 

• 

^ Shape memory alloys [SMA] are materials having many applications and attracting a 

^-H lot of interest due to their unique properties of shape memory and pseudoelasticity, which 

stem from both temperature-induced and stress-induced martensitic phase transition. The 
. ^ mechanical behaviour of such materials is rather complex and arises from a strong interaction 

^ between thermal and mechanical phenomena. In fact, when these materials undergo the 

^ martensitic phase transition, the increase in temperature due to localised self-heating/self- 

cooling has been experimentally found to be anything but negligible [H EJ [3l [4]. Hence, 
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the effects of heat transfer and of the heat dissipation towards the ambient are of primary 
relevance in the study of these materials [21 [5] and play a role in the rate-dependent behaviour 
of SMA [3, 6J. 

In modelling the mechanical behaviour of a SMA it is therefore important to account 
for mechanical and thermal aspects at the same time and to evaluate all the phenomena in 
their time evolution. 

In literature there exist many constitutive models for SMA [7J , which are derived follow- 
ing different approaches. Some authors propose models developed within plasticity frame- 
works, which make use of one or even more internal variables to describe the pseudo-elastic 
behaviour [HI [9l ^0\. Another group of constitutive models has been developed within a 
thermomechanic context; in his work, Tanaka [H] proposed a macroscale model comprising 
an internal variable to quantify the extent of the phase transition, a dissipation potential 
and an assumed transformation kinetic. Similarly, other models based on the same ther- 
modynamic background but featuring different kinetic laws have been formulated p!2l [T3]. 
and a complete heat equation which accounts for the contribution of latent heat and of 
the dissipation connected with the phase transition was added to the framework [HI [15] . 
Among others, Chang et al. [16] proposed a model comprising a strain-gradient elastic free 
energy and a chemical (transitional) free energy; Abeyaratne et al. [17] developed a model 
including a non-convex Helmholtz free energy function of the strain in order to describe the 
regions of stability of the different phases. In the context of using a non-convex free energy, 
it is not uncommon that micro-mechanical models are proposed [HI [19] ; MuUer and Seelecke 
fT9] developed a microscale model based on statistical physics comprising a bulk free energy 
function of the lattice shear deformation. 

Another approach to the description of the behaviour of SMA consists in applying the 
Ginzburg-Landau theory for phase transitions. At a microscopic scale, Falk [20j developed a 
Landau theory based on a shear strain order parameter. Levitas and Preston [21] proposed 
a single-grain model based on the decomposition of the strain in an elastic and a transforma- 
tional part, the latter being described by a pure order parameter which can assume different 
values to identify the different phases. In a similar fashion, Wang and Khachaturyan [22j 
proposed a model in which the strain field depends on several order parameters, whose 
evolution is described by several time-dependent Ginzburg-Landau [TDGL] equations. Nev- 
ertheless, attempts to handle the problem at a larger scale of observation have been made 
as well. Among others, Ahluwalia et al. [23j and Chen and Yang [24] proposed meso-scale 
models for the description of polycrystalline materials in which the order parameters ac- 
count for the different orientations of each single grain; in the same spirit, Brocca et al. 
[25] proposed a microplane model attempting to bridge the gap between micromechanics- 
based and macroscale models. Berti et al. [26j proposed a Ginzburg-Landau model which 
can be applied at the macroscale and encompasses mechanical as well as thermal effects by 
introducing the balance of linear momentum equation and the heat equation in a thermo- 
dynamically consistent framework; the order parameter is related to the extent of the phase 
transition between austenite and the martensite variants and its evolution is regulated by 
a TDGL equation. Owing to the presence of the heat equation, it is possible to describe 
the thermo-mechanical interactions which strongly influence the constitutive behaviour of 
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a SMA and to account for non-isothermal conditions. The model is presented both in a 
three-dimensional and a monodimensional setting; the reduction to this latter case is still 
meaningful because in most engineering applications SMA wires are employed. 

In this paper, we present a numerical study on the thermomechanical behaviour of a 
SMA, with a particular focus on the rate-dependent response and on the influence of ther- 
mal conduction and heat transfer in the mechanical behaviour. Starting from the models 
proposed in Berti et al. [26j and Daghia et al. [27j, we formulate a new free energy functional 
and give a proper expression for the relaxation parameter regulating the TDGL equation; 
this gives the model the properties required to reproduce the symmetrical behaviour between 
the austenite-to-martensite and the martensite-to-austenite phase transitions observed ex- 
perimentally. The ability of the model to quantitatively reproduce a variety of experimental 
evidences of a typical poly crystalline NiTi alloy is then demonstrated. 

The paper is organised as follows. In Section [2} the theoretical model is described, start- 
ing from a suitable Ginzburg-Landau free energy whose properties are outlined in Section 



2.1 Constitutive relations are given in Section [2^21 while the thermodynamic consistency 



of the model is shown in Section |2.3| and the complete differential system with appropriate 



boundary and initial conditions is provided in Section |2.4[ In Section |3] the Galerkin formu- 
lation of the differential problem, suitable for the flnite element implementation, is sketched. 
Section [4] is the main part of the paper, where the numerical results of a number of simulated 
tensile tests on a bar specimen under different conditions are reported. After describing, in 



Section 4.1, the effectiveness of the model in recovering the main experimental evidences in 



terms of stress-strain response, phase morphology and temperature evolution, in Section [42 
tensile tests performed at different values of the initial temperature are illustrated, while 
in Section |4.3| the stress-strain response under partial-loading conditions is depicted. The 
rate-dependent behaviour of the model is investigated in Section |44| with a particular focus 
on the effects of the strain-rate on the domain nucleation, the hysteresis cycle and the energy 
dissipation. The influence of heat transfer phenomena on the mechanical response of the 



specimen is examined in Section 4.5, The last two aspects are considered jointly in Section 



4.6 The paper ends drawing some conclusions in Section [5) 



2. Model 

The thermomechanical properties of SMA stem from the nature of the martensitic phase 
transition; here we will account for it by means of a Ginzburg-Landau approach [t20j, i21j . 
To this end we have to deflne a phase fleld, or order parameter, Lp] differently from the 
approach followed, among others, by Falk [20], the order parameter is not identifled with 
the uniaxial strain 6:11, but rather the phase fleld ip is used here as a macroscopic indicator 
of the phase (martensite or austenite) of the material at every point [211 ESI EHl EHl [30] . 
Strain is an independent fleld, though it is coupled with the order parameter. We assume 
that the martensite phase be present in only two relevant variants, M±, characterised by 
opposite uniaxial transformation strains with respect to the austenite phase [A], namely 
±St = ±8^)71 ® n, where n is a flxed direction. The Ginzburg-Landau potential, whose 



properties will be described in detail in Section 2.1, is expressly formulated to have (at most) 



three stable values of Lp at values {0, ±1}, corresponding to phases A, M±, respectively. We 
remark that, as a consequence of the first order character of the martensitic phase transition, 
intermediate values of ip do not represent any physically stable phase, but only transition 
layers between different stable or metastable physical phases. 

The model has to encompass evolution equations for the three basic fields: the order 
parameter (99), the strain (s) and the temperature {9). The evolution of the phase field will 
be described by the TDGL equation 

The strain will be determined through the combination of the balance of linear momen- 
tum equation and a proper constitutive relation between strain, stress and order parameter, 
while the heat equation regulating the temperature evolution will arise from the balance of 
energy equation. 

For the formulation of the model, we consider a bar sample of SMA occupying a bidi- 
mensional material domain ^. We indicate with x = (xi,X2) the material coordinates and 
with u{x) the displacement vector field. Relying on the small displacement approximation, 
we use the linearised strain tensor e = V'^u = |[Vu + (Vu)^]. 

2. 1 . Ginzburg- Landau potential 

Following Berti et al. [26j, we assume a Ginzburg-Landau free energy in the form 
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{Oa - Om) F{ip) +{9- ^l^stgniip)) G{ip)] , (2) 



where the temperature function 9 is given by 

e = max{(^ - Oa) , {do - 6 a)} , ^o < ^m < 0a, (3) 

being c the specific heat, 9ref a reference temperature for the specific heat, A the compliance 
modulus, a the stress tensor, k. the interface energy parameter regulating the interface 
width, and ^, 9a^ 9m ^ ^o material parameters related to the phase transition whose physical 
meaning will be explained in the following. 
The two potentials F and G are defined by 

F{ip) = l[l-cos(7r^)] + ^[l-cos(27r^)], (4) 

G{lp) = 1 — COS (nip) . (5) 

and are depicted in Fig. [l](a). The picture shows the shape of the potentials in the interval 
[—1, 1], which is the physically possible range of the order parameter. 













\ 








/ 


\ 








/ 


\ 








/ 


\ 








/ 


\ 








1 












\ 








1 


\ 








1 


\ 








1 


\ 








1 




\ 




1 


y^ 




\ 




1 


/ 




\\ 




1/ 




\\ 




1/ 






\\ 




// 













-1 -0.5 0.5 



\—f — g\ 

fa) 



\ 
\ 
\ 


\ 


y 

/ 
/ 

f 




\ 1 
. \ 

■••A 


/ 

/ y 
/y 




1 T^/'^^° 


~-'x^. 1 


^ 


/ 


\ ^. 


.y / 








\ 


/ 






\ 


/ 






\ 








\ 



5 = 0.5 5 = 5=-.5 s=-\ 

s=-\5 



(b) 



Figure 1: (a) Functions F{^) and G{^^) defining the free energy, (b) The potential W{^i^\s) in equation (|6| 
for different values of the parameter s. 



These potentials are indeed different from the polynomial ones usually employed in Lan- 
dau theory, in particular in Berti et al. [261 EOj. The choice we made here was suggested 
by a preliminary numerical investigation, which showed the drawbacks of using polynomial 
potentials. To clarify this issue, we first point out that the functions F and G are defined 
in a way that, in the interval [—1, 1], the function 



W{^-s) = F{^)^sG{ip) 



(6) 



has minima only at 99 = or at 99 = ±1 for any value of the parameter 5 G M. More precisely 
if 5 > 0, Vl^ has an (absolute) minimum in 99 = and two inflection points in 99 = ±1. For 
s < —\ there are two absolute minima at 99 = ±1 and a relative maximum at 99 = 0. If 
— 1 < 5 < there are three relative minima located at 99 = 0, ±1; in particular, for s = —1/2, 
the three minimal points have equal values. Similar properties can be achieved using sixth 
order polynomials. However, by using trigonometric functions it is possible to obtain a 
more symmetrical potential as it regards to the properties of the minima in the three points 
(f = 0, ±1. More precisely, using trigonometric potentials, the following property holds 



W{^ ± 1; -5 - 1) = W{^] s)-l- 2s. 



(7) 



In particular, W{(f ± 1; —1) = W{(p] 0) — 1 (see Fig. [T](b)); this means that, at the critical 
value s = —1^, the shape of the incipient maximum at 99 = is identical to that of the 
incipient maxima at 99 = ±1 for s = 0~ . As a consequence, the driving force for the phase 
transitions A -^ ±M and ±M -^ A are the same. This circumstance allows to obtain a 
symmetrical behaviour of the stress-strain curves in the loading stage with respect to the 
unloading stage which would not be possible to achieve with polynomial potentials. 

The structure of the minima of the free energy function ^ in the half-lines (f > and 
(/9 < is determined respectively by the parameters 
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s± 



Oa — Om 
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St : cr 



according to the following rules: 

- if 5± > 0, then Lp = ±1 is a relative maximum (unstable), 

- if 5± < 0, then (f = ±1 is a relative minimum (stable or metastable), 

- if 5+ < — 1 or 5_ < —1, then (/; = is a relative maximum or a inflection point 
(unstable), 

- if 5+ > — 1 and 5_ > — 1, then (/p = is a relative minimum (stable or metastable). 
Thus the equations 

give rise to piecewise straight lines in the 9 — a plane (see Fig. [2]). The temperature 9q has 




Figure 2: The equilibrium phase diagram resulting from local free energy minimisation. In the insets the 
potentials W{(p; s) in the range — 1 < (/? < 1 for the values of s corresponding to the selected region of the 
plane 6 — a are depicted. 

to be chosen lower than 9m ^ otherwise 9 > 9m at any temperature and the local minimum at 
(/9 = would never disappear (so the phase A would always be metastable, even at < 9m)' 
The three independent temperatures 9a^ 9m ^ ^o, as well as the energy I and the strain ^o 
are adjustable parameters for the transition. The ratio ^/^o, gives the slope of the straight 
lines in the a — 9 plane, while ^o is the transformation strain. 



2.2. Stress- strain constitutive relation and thermodynamic potentials 

After having described the Ginzburg Landau potential which appears in the phase field 
equation ([l]), we give the phase dependent stress-strain constitutive relation which comple- 
ments the balance of linear momentum equation. 

Following Levitas and Preston [21j and Berti et al. [26j, we assume the following consti- 
tutive relation: 



-^ = Ao- + ^sign{i^)G{^). 

OCT Z 



(8) 



The tensor St = e^n ® n describes a deformation in the n direction. From this relation, 
recalling that G{1) — G{0) = 2, it can be noted that St represents the strain difference 
between the phase ip = 1 and 99 = at constant stress {transformation strain). 

we noticed that I/sq corresponds to the slope of the transition lines in the 



In Section 



2.1 



stress-temperature phase space. We will see now that £ is also related to the latent heats of 
the transitions. To this end, we define the entropy function as 



"d9 



^ref 



-9 



£09 
239 



GM 



^ref 



-e--H{6-6o)G{^). 



(9) 



This definition follows from standard thermodynamics arguments, assuming that ip corre- 



sponds to the Gibbs free energy. In Section |2.3| we will show that this assumption and the 
subsequent expression of rj (as well as of the internal energy e) are in agreement with the 
second law of thermodynamics. The latent heat for the transition from an initial phase i to 
a final one / is 

Li^f = Otr {Vf - Vi) , (10) 

where 9tr > 9m is the transition temperature; given this definition, L > means that a 
positive amount of heat is absorbed by the system. For a = the transitions occur at 
temperatures 9a and 9m] in particular we have 



^A^M± 



-9 



M^ 



Lm±^a 



(11) 



As the Ginzburg-Landau potential ij) corresponds to the Gibbs free energy, the Helmholtz 
free energy is defined by 



ifj = ip + ere 



^ a2|^2|'^IV7l2- 






(12) 



We note that rj = —i/jq also holds. Finally, the internal energy is accordingly given by 



ZUj^Qf Z Z Z L 



{OA-OM)F{ip) + (e-ee')G{ip) 



(13) 



2.3. Thermodynamic consistency 

In this section, we show that the constitutive stress-strain relation and the expression of 
the thermodynamic potential are consistent with the second law of thermodynamics in the 
form of the Clausius-Duhem inequality 

^>-v^(3) + r. (14) 

where q is the heat flux and r the external heat supply. To this end, an appropriate form 
of the flrst law has to be stated: 

e = V; + Vi^-V-q + r, (15) 

where V^^ and P^ are the internal powers associated, respectively, to the balance of linear 
momentum equation 

pu = V-o- + pb, (16) 

and to the Ginzburg Landau equation ([l]). Exploiting the flrst law to eliminate the source 
term in the Clausius-Duhem inequality and recalling the deflnition of the Helmholtz free 
energy ip = e — 9r]^ the reduced inequality is obtained 

ij + r^e-v;-vi^ + ^-v9<o. (17) 

For the internal mechanical power we adopt the standard expression 

V^ = a:i. (18) 



The power balance associated to the Ginzburg-Landau equation is obtained by multiplying 
that equation by ip [3ll [321 [3^ : 

Tip'^ = — '0(^(/i^ + (/3V • {K\/ip) = — '0(^(/i — i^\^(f ' Vip + V • {K(p\/ip) (19) 

The internal power is then deflned by 

7^; = r(^2 + ^^^ + ^v^ . v^^ (20) 

while V^ = V • {n(p\/Lp) is the balancing external power. 

By substituting the power expressions and the Helmholtz free energy ip = ip -\- a : e in 
the reduced equality, we obtain the condition 

(iPe + 7i) e + (ijj^ - ey. & - Tip^ + ^ -Ve <{), (21) 

which has to hold for every process. This is satisfled by the given constitutive stress-strain 
relation e = ipa- and entropy relation rj = —ipo^ with the further conditions 

r > and q = -kV9, k > 0. (22) 



Using the expressions of the internal powers (18) and (20), we also write in an explicit way 
the balance of energy (15): 



g 



'%-^ 



ref 



00"G{(p)0 + ee'G{(p) -T^^ = kAe + r. (23) 
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2.4' Differential system and boundary conditions 

In this section we collect the equations of the model and state appropriate boundary and 
initial conditions for the problem. 

As it regards the balance of linear momentum equation, we will adopt the quasi-static 
approximation neglecting the inertial term pii. This corresponds to assume that the time 
scales of thermal phenomena and phase evolution are larger than the stress wave time scale. 
We also consider experimental conditions in which the external momentum source b and 
heat source r vanish. We obtain 

r((^)(^+| {eA-OM)F\^)+{e-'-fn®n\cTsign{^)\G\ip) - /^A(/; = 0, (24) 
c^ - i^ K^ - ^o)G(^)j 9-l9 H{9- 9o)G\^)^ - r^^ = kA9, 

with 

Act = V^u — — -n ® nsign{(f)G{(f), (25) 

r{(f) = Tmin + 2 i^rnax " Tmin) [l " COs(27r(/9)] . (26) 

The relaxation parameter r in the Ginzburg-Landau equation has been chosen on the 
basis of the experimental data available in literature [H El Hj. 

We now state the boundary and initial conditions accompanying the differential system. 
We consider a two-dimensional rectangular domain Q = [0, L] x [0,a]. With dQ^^dQ- 
and 5^1, 91^2 we denote two different partitions of the boundary dQ. At any point of the 
boundary, u denotes the outward unit normal. Typical boundary conditions are 

(27) 
(28) 
(29) 
(30) 
(31) 



\/(p{x,t) -U OQ 


= 0, 


U{x,t) 9Q_ 


= u{x), 


Vu{x,t) -u dn+ 


= tix), 


0{x,t) 9^1 


= e{x), 


V9ix,t)-iy9a, 


= gix). 



We remark that equation (27) represents the boundary 'insulation condition', which states 
that no long distance interactions are allowed between the domain and its exterior [34j; 
Dirichlet boundary conditions are not physically meaningful for the phase parameter. Mixed 
Dirichlet and Neumann boundary conditions are instead common for the temperature and 



the displacement fields. In the numerical tests, a different form for equation (31) will be 
used: 

kV9{x,t)'iy\9n, = h{9-9e.t). 

which is a convective heat transfer condition, representing a situation in which a bound- 
ary heat flux driven by the temperature difference with the external environment (via the 
exchange coefficient h) is prescribed. 
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As it regards the initial conditions, they have to be set only for the temperature and the 
phase variable: 

^(x,0)=eo(x), (^(x,0) = (^o(x), (32) 

as for the displacement field u we are using the quasi-static equation ([24|l). 



3. Galerkin formulation of the problem 

In this section, a formulation of the governing equations is presented which will be used 
in computing numerical examples. 

To formulate the finite element problem, let T = {E} be a triangulation of the domain 
Q into finite element cells. We will work with the usual Lagrange finite element basis with 
the definition of the following piecewise polynomial spaces: 

Wk = {v^ G H' (n) , v^\^ G Pk {E)\/E G T} , 

Vk = {v^ eWk:v^ = on dft} , 

where Pk denotes the space of Lagrange polynomials. The numerical formulation that we 
propose involves the solution of the following problem: given the data at time t^, find 

(u^, (f^, e^) G {Wkf ^Wk^^Wk (where fc > 0) at time t^+i such that 



L^h 



L^h 



^ Vw'^ • i (Vu^:, - ^sign{^l)G{^l)) dx = Vw'^ e [V^f 



'^ 2 



{Oa - Om) F'iipt) +{e- ^^szgniipt) ) G'{^t 



dx + 



(33) 
(34) 



/ 



'a 



h^^n+l 






LqH — 



z 



9, 



ref 



-A Kot - Oo)G{^t) 



nh nh 



. ^n+1 - ^n 



At 



dx 



dx 



(35) 



Vz^ • k {ve'')^dx + 






where (.)m = (1 - /5)(-)n+i + /5(-)n. 

Ly/i, L^h and L^h represent the weak forms of the three governing equations: the balance 
of linear momentum, the TDGL equation and heat equation, respectively. The forms are 
cast, following a standard process, using the Crank-Nicolson method for the time derivatives; 
as the functional L^^, L^h and Lqh are nonlinear in (i6^, Lp^^O^)^ a modified Newton method 
is used to solve the nonlinear problem at each time step. In the numerical simulations, 
linear Lagrange functions are used for the phase and the temperature fields, while quadratic 
Lagrange functions are used for the displacement field. 

The computer code used to perform all the simulations has been generated automatically 
from a high-level implementation by using a number of tools from the FEniCS Project |35j. 
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Figure 3: Imposed displacement at the end of the bar as a function of time. 

4. Numerical simulations 

A number of numerical simulations were performed with the aim to describe the effec- 
tiveness of the model in capturing the behaviour of a SMA not only from a qualitative, but 
also from a quantitative point of view. The results obtained are presented in this section 
and their consistency with the experimental evidences available in literature is highlighted. 

To reproduce the results of a tensile test conducted on a SMA wire or on a dog-bone 
specimen, a displacement-controlled uniaxial test was simulated considering a domain of 
dimensions L = 30 mm and a = 0.08 mm. The n direction appearing in the transition 
tensor St is assumed to be along the major domain edge L. The imposed displacement 
varied linearly from ti = to a maximum value u = Umax'^ (loading stage) and then back 
to 'U = (unloading stage), as depicted in Fig. [3} In all cases, the initial temperature is set 
equal to the external temperature {O^xt = ^o) 

To explore the potential of the model in capturing the experimental evidences, the values 
of the parameters used in the simulations were chosen according to the data reported in the 
paper by He and Sun |4j regarding a commercial polycrystalline NiTi sheet. The parameters 
for which it was not possible to directly find a value in the aforementioned paper were 
chosen in order to reproduce in a satisfactory way the response observed by Zhang et al. [3j, 
especially in terms of transition stress, propagation stress, number of nucleating domains 
and maximum increase in temperature during the test. 

The values of the parameters are: A = 2.5 x 10"^^ Pa"^; ^o = 0.058; / = 10^ Pa/K; 
Oa = 288.5 K; 9m = 273.0 K; ^o = 212.7 K; /^ = 0.15 N; Tmax = 8 x 10^ Pa s; Trmn = 
1.6 X 10^ Pa s; c = 3.2 x 10^ Pa/K; 9ref = 296 K; A; = 18 W/(m K). 

The values of the heat transfer coefficient /i, of the initial temperature Oq and of the 
imposed strain rate e varied depending on the test. 

^.7. General features of a tensile test 

The nominal stress-strain diagram resulting form the simulation of a tensile test is re- 
ported in Fig. |4](a), with the related microstructure evolution depicted in Fig. [4](b). Sim- 
ulations are performed at a nominal strain-rate e = 3.3 x 10~^ s~^ and with a value of the 
heat transfer coefficient h = 10 W/(m K) at an initial temperature 9q = 296 K. 
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Figure 4: Tensile test at nominal strain rate e = 3.3 x 10~^ s~^, initial temperature ^o = 296 K and heat 
transfer coefficient h = 10 W/(m K): (a) nominal stress-strain diagram; (b) phase morphology at different 
stages. 



The specimen is initially in an unstressed configuration and in the austenitic phase. When 
the test starts, the specimen is stretched and austenite deforms elastically; the deformation of 
the domain is therefore homogeneous. This behaviour is maintained until the domain reaches 
the point A in Fig. |4| which is the configuration at which austenite becomes unstable and 
martensite starts nucleating. At this strain rate, three martensite spots arise; the nucleation 
stage is completed at point B in Fig. [4| when three martensite bands are clearly visible. 
According to Chang et al. [16], we remark that within this context the term nucleation 
refers to a macroscopic process, i.e. the origination of the transformation fronts associated 
with the macroscopic localized transformation, rather than nucleation at microscopic scale 
as it is properly meant in material science. The main results of the nucleation process are a 
localised deformation of the martensitic bands which causes the stress to drop by 20 MPa 
and the formation of a number of transition fronts. 
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After originating, the fronts start propagating along the domain; this stage corresponds 
to the plateau of the stress-strain curve. Interestingly, a plot of the bulk transitional free 
energy function for the stress and temperature occurring at the transition fronts during 
phase propagation reveals that the transformation, once it has started, is able to proceed 
even if there is an energy barrier between the local minimum and the global minimum. 
Owing to the presence of the gradient term in the free energy, the points of the domain in 
the local minimum (A phase during loading, M+ phase during unloading), can overcome 
the energy barrier and evolve towards the configuration of absolute minimum (M+ phase 
during loading, A phase during unloading). The progress of the transition in the plateau 
regime is regulated by the interplay between thermal and mechanical phenomena. The heat 
released during the phase change at the transition front leads to an increase in temper- 
ature; consequently, the energy barrier increases and the transition is inhibited. On the 
other hand, an increase in stress causes the energy barrier to decrease and the transition is 
favoured. The interaction between these two opposite trends sets the slope of the plateau 
in the stress-strain diagram, which may change considerably under diflFerent strain-rates or 
in different environmental conditions, as will be shown in Sections |4.4| and [45} During the 
phase propagation, front-merging phenomena may occasionally occur when two transition 
fronts merge together (points D and E in Fig. [4]) or when one front reaches the end of the 
specimen (point C in Fig. [4l the left transition front reaching the left end of the domain). 
These events are refiected in a sudden drop and rise of the nominal stress. 

At point E in Fig. [4] the specimen is fully transformed into one variant of martensite 
and retrieves an elastic behaviour (from E to F) in which the deformation of the domain is 
homogeneous, corresponding to an elastic martensite deformation. 

The M+ -^ A transition starts at point F in Fig. |4} due to the choice of the free 
energy function and of the relaxation parameter of the TDGL equation (see Section 2.1), 
it has analogous features to the A -^ M+ phase transition (note that the stress rise at the 
beginning of the unloading plateau is of the same magnitude of the corresponding stress 
drop in the loading stage and the number of product-phase nuclei is the same). The test 
ends with a residual deformation due to some residual martensite in the specimen. 

The A -^ M+ transition is exothermic, hence the specimen experiences self-heating in 
the neighbourhood of the transition fronts due to both latent heat release and dissipation; 
furthermore, there is an overall increase in temperature due to thermal conduction within 
the domain. Opposite considerations apply to the case of the M+ -^ A transformation. 
The temperature distribution along the specimen is shown in Fig. [5] for some of the points 
marked on the stress-strain diagram in Fig. |4} The curve marked with A refers to the onset 
of the transition; until this event, the domain has deformed elastically, hence there is no 
increase in temperature with respect to the initial condition. After the nucleation, three 
temperature peaks rise, as depicted in curve (7; the temperature profile exhibits then as 
many fronts as the phase transition fronts, because of the heat released during the phase 
change. When two propagating fronts merge, the released heat becomes more significant 
and the local increase in temperature is more pronounced; this is the situation depicted in 
curve E in Fig. [5| where the maximum temperature of the test is experienced. The values 
are quantitatively in good agreement with the experimental observations [4J. The opposite 
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behaviour is observed in the M+ -^ A transformation (curves G, H and L in Fig. [5]) during 
which the temperature decreases even below the initial temperature as a consequence of the 
latent heat absorbed in the transition from martensite to austenite. 
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Figure 5: Tensile test at nominal strain rate e = 3.3 x 10 ^ s ^, initial temperature ^o = 296 K and heat 
transfer coefficient /i = 10 W/(m K): temperature distribution along the specimen (letters refer to Fig. [4|. 



J^.2. Tensile tests at different values of the initial temperature 

The nominal stress-strain diagrams depicted in Fig. [6] show the results of tensile tests 
simulated with different values of the initial temperature, which ranged from Oq = 296 K to 
^0 = 373 K. 

According to the phase diagram in Fig. [2| a higher value of the temperature implies 
a higher stress required for the macroscopic domain nucleation. As a result, Fig. |6] shows 
that as the initial temperature is higher, the hysteresis cycle appears shifted towards higher 
values of stress and the elastic behaviour is retrieved after the hysteresis loop is completed. 
Note that the area of the hysteresis loops, as well as the number of product-phase nuclei 
arising at the onset of the transition are independent of the value of the initial temperature. 



4^3. Tensile tests at different values of the maximum deformation 

To evaluate the behaviour of the specimen under partial-loading conditions, tensile tests 
are performed at different values of the maximum imposed nominal strain, which ranged 
from e = 0.07 to e = 0.027. Fig. [7] shows the resulting stress-strain diagrams and the 
morphology of the specimen at interesting stages of the tests. 

It can be observed that if the unloading begins before any domain merging, the reverse 
transformation consists in a propagation of the already formed transition fronts, rather 
than involving a re-nucleation process (cases A, 5, C and D of Fig. [7|. If the loading 
is otherwise interrupted after there has been a merging of two fronts, there can be re- 
nucleation. Considering case E of Fig. [7[ at the end of the loading stage there is only 
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Figure 6: Tensile tests at nominal strain rate e = 3.3xl0^s^ and heat transfer coefficient h 
10 W/(m K): nominal stress-strain diagram for different values of the initial temperature. 



one austenite domain left; the reverse transition starts therefore with the nucleation of one 
austenite domain (£"') and the transformation proceeds with two propagating domains. 

Case F, instead, depicts the situation in which the specimen is fully transformed into 
martensite upon loading; as a consequence, the reverse transition must start with the nucle- 
ation of austenite domains. As can be seen in Fig. [7|(b), in this case three nuclei originate 
and consequently the rise in stress due to the nucleation is higher than in case E. 

In much of the data available in literature regarding SMA, the specimen is fully trans- 
formed into martensite only in the gauge section; two austenitic domains remain at the two 
ends of the specimen, and the reverse transformation starts with the propagation of the two 
existing fronts without re- nucleation [3j . The mechanical response of the specimen suffers 
therefore from an experimental artifact and the onset of the reverse transformation does not 
exhibit any stress rise due to nucleation. This condition is replicated in cases form Aio D 
of Fig. [7| To overcome this issue and ensure product-phase nucleation in both loading and 
unloading stages, ladicola and Shaw [36j make use of a thermoelectric device to adjust the 
temperature at the clamped ends of the specimen in order to ensure the transition to start 
within the gauge length and obtain a stress-strain response similar to that described in case 
F of Fig. [7| 
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Figure 7: Tensile tests at nominal strain rate e = 3.3 x 10~^ s~^, initial temperature ^o = 296 K and 
heat transfer coefficient h = 10 W/(m K) for different values of the maximum deformation: (a) nominal 
stress- strain diagram; (b) phase morphology at different stages. 



4^4' Influence of the strain-rate on the behaviour of the model 

To explore the ability of the proposed model to reproduce the rate-dependent behaviour 
experimentally observed in SMA j3], tensile tests at different values of the nominal imposed 
strain rate were performed; the value of the strain rate ranged form e = 3.3 x 10~^ s~^ to 
e = 3.3 X 10~^ s~^, while the initial temperature and the maximum imposed displacement 
were the same for all the test-cases (^o = 296 K, e = 0.07). 

As the strain-rate becomes higher, the stress at which nucleation starts increases (Fig. 
[sFa)); as a consequence, the driving force for the phase transition is greater and more 
product-phase nuclei will originate (Fig. |8](b)). Fig. [9](a) depicts the maximum number of 
product phase nuclei as a function of the nominal imposed strain rate; the trend resulting 
from the simulations is exponential, in agreement with the experimental observation of He 
and Sun [4] and the numerical simulations of ladicola and Shaw j37j . 

Fig. |8](a) also shows that the discontinuities in the stress-strain curve induced by nucle- 
ation and front merging tend to vanish at high strain rates, because the loading time scale 
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Figure 8: Tensile tests at different values of nominal strain rate, initial temperature 6>o = 296 K and heat 
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10 W/(m K): (a) nominal stress-strain diagram; (b) phase morphology after M^ 



gets closer to the time scale of the phase transition phenomena. 

One of the main effects of the strain-rate on the mechanical behaviour of the specimen 
is the change in the slope of the stress plateau during the loading stage, which results from 
an interplay between transitional, thermal and mechanical phenomena. As the strain rate 
becomes higher, the heat released due to the phase transition increases, but there is less time 
for such heat to be transferred to the exterior, because the duration of the loading stage 
is shorter; as a consequence, at a given strain the temperature increases. In order for the 
transition to proceed, the stress must therefore increase, which makes the slope of the stress 
plateau increase. Moreover, if the stress at the end of the loading stage is higher, the elastic 
deformation of the specimen is greater, hence less deformation from the phase transition is 
required to achieve the maximum imposed deformation and more untransformed austenite 
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Figure 9: (a) maximum increase in temperature within the specimen and number of product-phase nuclei 
as a function of the test strain rate, (b) Energy dissipation as a function of the test strain-rate. 



remains. Furthermore, upon unloading the stress rise described in Section 4J_ is absent 
at higher stress rate; this effect is due to the fact that at the fixed maximum strain there 
remains more untransformed austenite, hence no re-nucleation occurs. 

As a direct evidence of the scarcity of time for the heat generated within the specimen 
to be transferred outside it, the maximum temperature reached at the end of the loading 
stage is higher for higher strain rates. The values of the maximum temperature increase 
within the specimen are reported in Fig. |9} in which it can be noted that higher strain-rates 
correspond to greater temperature increases. 

The strain-rate has a strong impact also on the damping capacity of the specimen; a 
measure of this property is the area of the hysteresis cycle in the stress-strain diagram, 
which represents the mechanical energy dissipated during the test. Fig. |9](b) shows that the 
dissipated energy changes non-monotonically with the strain rate. As observed by Zhang et 
al. [3j, the non-monotone trend of the dissipated energy is due to the fact that during the 
A -^ M+ transformation the slope of the stress plateau increases monotonically with the 
strain rate, hence the transition occurs for higher values of the stress; on the other hand, 
the M+ -^ A transition has a non-monotone trend [3j . 

4^5. Influence of the heat transfer conditions on the behaviour of the model 

The rate-dependent behaviour of SMA stems from an interaction between transitional, 
thermal and mechanical phenomena; hence, changing the conditions under which the gen- 
erated heat is exchanged with the exterior affects the mechanical response of the material 
[5j. To highlight this point, tensile tests simulating different heat exchange conditions were 
performed; this was achieved by changing the value of the heat transfer coefficient, which 
ranged form h = 2 W/(m K) to h = 50 W/(m K). For all the test-cases, the initial 
temperature and the imposed strain rate were the same (^o = 296 K, e = 3.3 x 10~^ s~^). 



The nominal stress-strain diagram for the different tests is reported in Fig. [T0| A low 
value of the heat transfer coefficient causes an increase in the time needed for the generated 
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Figure 10: Tensile tests at nominal strain rate e = 3.3 x 10 ^ s ^ and initial temperature ^o = 296 K: 
nominal stress-strain diagram for different values of the heat transfer coefficient. 



heat to be transferred to the exterior. If the value of the coefficient is very low, the test 
proceeds in adiabatic condition, as the time scale of the heat transfer phenomena becomes 
much larger than the test time scale. As a result, in the loading stage the temperature at 
the transition front increases considerably and a higher stress is required for the transition 
to proceed. On the contrary, if the heat transfer coefficient is high, the generated heat is 
transferred to the exterior in short time and the test proceeds nearly isothermally; under 
this condition, the temperature does not contribute to increase the energy barrier and the 
transition can therefore proceed at a low value of the stress. As a result, the stress plateau 
is horizontal and the difference between the nucleation stress and the propagation stress is 
relevant. 

J^.6. Concurrent effect of ambient condition and strain-rate on the behaviour of the model 

So far, the effects of the strain-rate and of the heat transfer on the mechanical behaviour 
of SMA have been investigated disjointly; in this section, we discuss the concurrent inffuence 
of these two aspects on the damping capacity and on the maximum temperature experienced 
during the tests. 

Fig. [TT](a) depicts the maximum increase in temperature during the simulated tensile test 
for different values of the heat transfer coefficient and of the imposed strain rate. Both a high 
value of the strain rate and a low value of the heat transfer coefficient cause the maximum 
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Figure 11: (a) maximum increase in temperature within the specimen and (b) energy dissipation as a 
function of the test strain-rate for different values of the heat transfer coefficient. 



temperature inside the specimen to increase, as a smaller amount of the generated heat 
can be transferred to the exterior. Note also that at high strain rates the temperature 
increase within the specimen tends to become insensitive to the heat transfer coefficient, as 
the adiabatic condition is approached regardless heat transfer phenomena. 

The non-monotone behaviour of the damping capacity for different values of the strain 



rate was analysed in Section |4.4t the same test was performed for different values of the 
heat transfer coefficient {h = 5 W/(m K), h = 20 W/(m K)) and the results obtained are 
shown in Fig. [TT](b). It can be observed that a high value of the heat exchange coefficient 
causes the peak of the dissipated energy to be located at higher values of the strain rate, 
owing to the fact that heat transfer towards the exterior is favoured and the adiabatic 
condition is shifted towards higher values of the strain rate, thus. Note that, for the range 
of values of interest in this study, there is no significant effect of the heat transfer coefficient 
on the maximum value of the dissipated energy. 

The stress-strain diagrams referring to the tests performed at different strain rates for the 
values of the heat transfer coefficient h = 5 W/(m K) and h = 20 W/(m K) are reported 
in Figs. [12] and [131 respectively, together with the maximum number of martensite nuclei 



originated during the loading stage. The behaviour is the same described in Section [44] for 
the case h= 10 W/(m^ K). 



5. Conclusions 

A Ginzburg-Landau model that reproduces the most relevant macroscopic features of 
the behaviour of a shape memory alloy has been considered. The model encompasses a 
time-dependent Ginzburg-Landau equation for the evolution of the phase order parameter, 
the balance of linear momentum and the heat equation, in order to account for the interplay 
between the transitional, mechanical and thermal aspects which strongly influences the 
behaviour of SMA. Starting from the theoretical model described in Berti et al. |26], a new 
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5 W/(m K): (a) nominal stress-strain diagram; (b) phase morphology after M^ 



free energy based on trigonometric functions has been formulated and a suitable form of the 
relaxation parameter of the TDGL equation has been proposed. The resulting properties of 
the model ensure the driving force for the phase transition from austenite to martensite to be 
equal to the driving force for the transition from martensite to austenite and, consequently, 
provide a symmetrical behaviour of the model upon loading and unloading. 

Numerical simulations have been performed to investigate the mechanical response of 
the model; tensile tests on a bar specimen have been simulated under different conditions 
by changing the initial temperature, the maximum imposed strain, the nominal imposed 
strain-rate and the environmental conditions. The results obtained illustrate that the model 
is able to reproduce the main experimental evidences reported in literature [H EJ [3l HJ [5] . The 
simulated tensile tests show the presence of a stress drop during the phase nucleation, which 
may be absent in the unloading stage if the loading is interrupted before the completion of the 
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Figure 13: Tensile tests at different values of nominal strain rate, initial temperature ^o = 296 K and heat 
transfer coefficient h = 20 W/(m K): (a) nominal stress-strain diagram; (b) phase morphology after M+ 
nucleation. 

phase transition; temperature profiles consistent with the observed physical phenomena have 
also been described. The effect of the strain-rate on the mechanical response of the model 
have been investigated; its impact on the number of nucleating domains, the nucleation 
stress, the stress relaxation during nucleation and on the slope of the stress plateau (strain- 
hardening) have been recovered. The role of the environmental conditions in the mechanical 
response have also been examined, and the concurrent infiuence of heat transfer conditions 
and strain-rate on the damping capacity and on the overall behaviour of the model have 
been described. 
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